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Abstract. The microturbulent approximation of turbulent motions is widely used in radiative transfer calculations. 
Mainly motivated by its simple computational application it is probably in many cases an oversimplified treatment 
of the dynamical processes involved. This aspect is in particular important in the analysis of maser lines, since 
the strong amplification of radiation leads to a sensitive dependence of the radiation field on the overall velocity 
\& ' structure. To demonstrate the influence of large scale motions on the formation of maser lines we present a simple 

stochastic model which takes velocity correlations into account. For a quantitative analysis of correlation effects, 
we generate in a Monte Carlo simulation individual realizations of a turbulent velocity field along a line of sight. 
Depending on the size of the velocity correlation length we find huge deviations between the resulting random 
profiles in respect of line shape, intensity and position of single spectral components. Finally, we simulate the 
O !' emission of extended maser sources. A qualitative comparison with observed masers associated with star forming 

regions shows that our model can reproduce the observed general spectral characteristics. We also investigate 
shortly, how the spectra are effected when a systematic velocity field (simulating expansion) is superposed on the 
fluctuations. Our results convincingly demonstrate that hydrodynamical motions are of great importance for the 
, understanding of cosmic masers. 
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F~ i 1. Introduction the hydrodynamical velocity changes substantially, plays 

an important role. In the usual microturbulent analysis 
Maser lines of various molecular species are observed in large gcale motions are neg lected. As a consequence, each 
star forming regions and in the envelopes of evolved late spectral component will generally be attributed to a sep- 
type stars. In both cases it is well known that turbulent arate source and CQuld in principle be formed under dif _ 
motions play an important role in the gas hydrodynam- ferent physical conditions. The aim of the present inves- 
ics. Their spectra show large velocity differences (up to tigation is to emphasize the necessity to consider hydro- 
tens of km/s) and are variable on observable time scales. dynamica l motions in more detail. We will show how cor- 
VLBI-observations reveal that maser regions often consist related veloc ity fields affect the line forming process and 
of several isolated sources with different radial velocities that om itting this important aspect can lead to erroneous 
(see, e.g., Reid et al. |1980). Due to the Doppler effect, interpretations of observed maser spectra. 

the radiation field depends on the local velocity struc- • i n , , lr , , , n hhA i t> 

_, , ,. , . . „ , , , ... Iravmg and collaborators (Gail et al. 1974, Iravmg 

ture. Ihis coupling to the velocity held does sensitively n t\h\± i. j i j ^ n ± r 

_ , , . ° , .. . ., 19751 developed a theory allowing to account for corre- 

aftect the formation of maser lines since m the unsatu- , , ; , , , , ,. i . r rn . ,. ,. 

. lated turbulent motions m trie transfer ot line radiation, 
rated regime the radiation is exponentially amplified. To 

the extent that stochastic hydrodynamical velocities are i TT 7~. 7, 77, • , , , , ■ 

J J Here and m the following we use the term turbulence m the 

involved, the correlation length, i.e., the length over which spectroscopic sense (which is wider than in hydrodynamics) 

referring to any stochastic motion leading to line broadening. 

Send offprint requests to: R. Boger The terms micro-, meso-, and macroturbulence refer to the 
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In their formalism, the hydrodynamic velocity v along a 
line of sight is described by a Markovian process defined 
by a Gaussian distribution function and an exponential 
correlation function. Due to the stochastic modeling of 
the velocity field, the intensity also has to be considered 
as a random variable. Traving showed that in accord with 
these assumptions and together with the usual equation of 
transfer, a Fokker-Planck equation can be derived, which 
describes the probability W{I V , v, s) of finding at point s, 
the intensity l v and the velocity v. Of particular interest is 
the expectation value of the intensity (J„) since it reflects 
as a time or spatial average the mean physical properties 
of the system. Gail et al. I|1975(l examined the model of an 
unsaturated maser applying the theory described above. 

While this approach gives a first qualitative insight 
into the effects introduced by a turbulent velocity field, 
the question arises to which extent the calculated expec- 
tation value of the intensity may be directly compared 
with observed spectra. The answer depends in particular 
on two issues. The first is the question of how many sta- 
tistically independent lines of sight contribute to a given 
observation. This obviously is a function of the spatial 
resolution. The second point is the question of how many 
spectra corresponding to individual lines of sight have to 
be averaged in order to obtain the expectation value with 
a certain accuracy. As will be shown later, in the case of 
unsaturated masers this number is very large. In order to 
investigate these questions we set up a Monte Carlo (MC) 
simulation scheme to generate individual realizations of 
the random velocity field. The calculated line profiles re- 
flect the stochastic nature of the turbulent motion along 
a given line of sight and allow a quantitative analysis of 
correlation effects. 



2. Model assumptions and the Monte Carlo 
method 

In the case of an unsaturated maser the occupation num- 
bers of the maser levels are essentially not affected by 
the intensity in the maser line. Consequently, the inten- 
sity variation along any line of sight depends just on the 
local physical conditions. In order to focus on the effect 
of correlated velocity fields, we assume a plane parallel 
slab with a homogeneous density and temperature dis- 
tribution. Additionally, we take the velocity field to be 
quasi-static, presuming that the dynamical time scales are 
large in comparison with the duration of individual obser- 
vations. The variation of the turbulent velocity along each 
line of sight is assumed to follow a Markov process and the 
projected velocity v(s) along a line of sight is treated as a 
random function of s. The one-point distribution function 



way in which the turbulent motions affect the line broadening 
which depends on the ratio of the mean free path of a photon 
to the correlation length of the velocity field. 
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Fig. 1. (a) Individual realizations of the one-dimensional 
velocity field, and (b) the corresponding distribution func- 
tion P (v) for the microturbulent (I = 0, dashed his- 
togram) and a mesoturbulent (l/L = 0.1, solid histogram) 
model 

of v is considered to be Gaussian with the rms turbulent 
velocity at 

„,2 



W(v) = 



1 



2-KGt 



exp 



2a\ 



(1) 



The conditional probability of finding at point s + As the 
velocity v + Av is 



P(v + Av,s + As | v,s) = 

1 exp J [Av + v(l-f)f 

^ (1 - P) P l 2a? (I-/*) 

with an exponential type correlation function 



/ (As) = exp 



I 



(2) 



(3) 



The correlation length I defines the length scale of the 
stochastic velocity variation. From Eq. J2J it follows that 
the expectation value of Av is 



(Av) 



1 — exp 



\As\ 
I 



and the variance is given by 



Var (Av) = of 



1 — exp 



2|As 



(4) 



(5) 
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The assumption that the variation of the velocity field 
along the line of sight may be described by a Markov pro- 
cess (implying an exponential correlation function) is the 
simplest way to introduce correlation effects into the the- 
ory of radiative transfer (see e.g. Gail et al. 1974, 1980). 
This assumption implies that only the one-point and the 
two-point distribution functions are important. All other 
correlation functions require the knowledge of higher order 
multi-point distribution functions. For later comparison 
with the work of others, we note that the (spatial) power 
spectrum of the velocity field is the Fourier transform of 
the correlation function. This means that (3) implies a 
one-dimensional power spectrum 



I 



1 + l 2 k 2 



(6) 



which for l 2 k 2 ^> 1 is well approximated by a power law, 
P(k) ~ k~ 2 . If we consider (3) as correlation function in 
three dimensions, i.e. if we consider s as a vector, we find 
for the three-dimensional power spectrum 



I 



(i + i 2 k 2 y 



(7) 



The choice of the correlation function (3) makes it rel- 
atively easy to calculate individual realizations of the ve- 
locity field (see below), as well as the expectation value 
of the line profile using the generalized radiative transfer 
equation (Gail et al. 1974, 1975). 

In accord with Eq. @ the distribution of the random 
projected velocities can be generated in a Monte Carlo 
procedure (Levshakov et al. H997fl : 



v( S + As)=t;Ja 2 (l~f 2 ) + fv( S ) 



(8) 



where £ is a normally distributed random number with 
(£) = and Var(£) = 1. The velocity at the edge of the 
cloud is given by v = cr t £(s = 0). 

The line forming process depends sensitively on the ra- 
tio <J t /v t h (see, e.g., Levshakov & Kegel ITOMjl . In the case 
"t/^th < 1) the line broadening is thermally dominated 
and the turbulent velocity correlation is insignificant. 
Contrary, if a t /vth > 1 correlation effects become pro- 
nounced. An inappropriately chosen step size can smear 
out the correlated structure of a random velocity field. 
This is avoided by satisfying the inequality 



Var (Av) < v 2 h . 



(9) 



From JSJ it follows that for at > t>th this is equivalent to 



I 



< In 



1 - (vth/vtY 



-1/2 



(10) 



To illustrate the influence of a finite correlation length, 
Fig. shows one realization of a completely uncorrclated 
(microturbulent) velocity field (I — 0, dashed line) and 
one realization of a mcsoturbulcnt velocity field (/ > 0, 
solid line). The lower panel represents the corresponding 
projected velocity distributions P (v) . As expected, the ve- 
locity in the limiting case is normally distributed, whereas 
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Fig. 2. Upper panels: Line profiles for two individual real- 
izations of the velocity field (solid curve) for the parameter 

set l/L = 0.1,<7t/i>th = 10> an< ^ T ° = — HO- For compar- 
ison also the spectrum of the intensity expectation value 
is shown (dashed curve). Lower panels: The optical depth 
of the random spectra 

the correlated velocity distribution deviates significantly 
from a Gaussian. This is due to the fact that for a finite 
ratio 1 1 L the statistical base represents an incomplete sta- 
tistical sample. Averaging over many realizations (either 
temporal or spatial) would yield an complete ensemble 
with a unique Gaussian distribution function. For sim- 
plicity, we neglect the dependence of the occupation num- 
bers of the involved maser levels on the radiation in the 
maser line. If furthermore the source function S u is as- 
sumed to be constant over the whole region, the equation 
of radiative transfer can be solved analytically. Neglecting 
the background radiation yields 

I v (s = 0) = S v {1 - exp \-r v (L)}} , (11) 
with the optical depth 



t v (s) 



(12) 



and the monochromatic absorption coefficient 

K v = kq<$(Av,v) , (13) 

$ (Av, v) being the local profile function. It should be 
noted that for population inversion the source function 
and the absorption coefficient are both negative. If one 
considers thermal Doppler broadening as the dominant 
effect, the local profile function can be written as 



$ (Av,v) 



1 



y/irAvn 



exp ■ 



(v - vo) 



Az^d 



(14) 
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Fig. 3. Convergence of the ensemble average of N random 
profiles (solid line) to its expectation value (dashed line). 
The parameter configuration is the same as in Fig. [2] 

where Aj/d is the thermal Doppler width and Vth the ther- 
mal velocity. Equally, the projected random velocity dis- 
tribution along a line of sight P (v) can be convolved with 
the thermal profile function to derive the optical depth 

oo 

Tv (L) = k L I P (v) $ (Ai/, v) dv . (15) 



3. Results 

Our model calculations mainly depend on three param- 
eters: the optical thickness of the slab To, the turbulent 
velocity cr t and the correlation length I. Due to the ran- 
dom nature of the velocity field, the optical thickness can 
not be defined uniquely. Therefore, it is characterized by 
the optical depth in the line center for vanishing turbulent 
motions, a-t = 0, which is then given by kqL. 

Figure demonstrates the influence of velocity cor- 
relations on the line forming process. The upper panels 
show the line profiles for two individual realizations of the 
velocity field along a given line of sight calculated with 
the same statistical parameters (solid curves). For com- 
parison, also the expectation values of the intensity are 
shown (dashed curves). The lower panels give the corre- 
sponding optical depths t v . In our simple model (constant 
density and constant pumping efficiency), the latter reflect 
directly the complex structure of the underlying velocity 
field (see Eq. m]). Due to the exponential amplification 
the intensity profiles appear to have a much simpler struc- 
ture than the T^-profiles. This is in particular true if a 
lower intensity cut-off is introduced. Therefore do high in- 



tensity spectra often consist of a single line only, which 
may exhibit a pronounced line shift. 

As is obvious from Fig. (see also Table the line 
profiles calculated for individual realizations of the veloc- 
ity field can differ substantially from each other and from 
the expectation value. The degree to which an individual 
realization can deviate from the expectation value can sta- 
tistically be described by the convergence behavior of the 
ensemble average. We find that in general the mean value 
approaches the expectation value very slowly. For the cho- 
sen parameter configuration, Fig. [3] shows a sequence of 
averages over N realizations. Even the superposition of 
5 x 10 5 individual spectra is not sufficient to match the 
spectrum of the intensity expectation value. For a smaller 
correlation length or larger turbulent motions the conver- 
gence is improved but still a substantial number of pro- 
files is required. This finding is contrary to the case of 
pure absorption. In the context of quasar absorption lines 
Levshakov et al. (1997) showed that the average of about 
100 realizations reproduces the expectation value of the 
intensity within a one percent error margin. The radia- 
tion of an extended source reflects the average influence 
of the dynamical motions on the radiation field. From the 
above finding it can be concluded that in the case of an 
unsaturated maser the expectation value of the intensity 
will in general not be in good agreement with the mea- 
surements if a finite correlation length is involved. Even 
if the problem is solved seli-consistently, where, due to 
saturation effects, the deviations from the mean intensity 
will be less pronounced, it can be expected that the mean 
intensity will not match the observations. Despite of this 
conclusion, the method by Traving is a valuable tool to ex- 
amine the general influence of the stochastic parameters 
which govern the dynamical motions. It allows a very fast 
calculation of expectation values and distribution func- 
tions (Levshakov et al. H997[) . In this respect, the method 
is preferable to a MC simulation which requires the gen- 
eration of a very large number of realizations and is thus 
costly in computer time. Starting from the same model 
assumptions, Gail et al. (1975) examined the mean prop- 
erties of the radiation field. Their most important finding 
is that the mean intensity increases with increasing veloc- 
ity correlation. This corresponds to our result that with 
a larger correlation length the probability of having large 
peak intensities grows (see Table In the present pa- 
per we want to focus on a different aspect of the problem 
for which the MC technique allows a more quantitative 
comparison of observations and model calculations. In the 
framework of our simple model, it is reasonable to ap- 
proximate the radiation of an extended maser region by 
the superposition of several spectra corresponding to sta- 
tistically independent lines of sight, i.e. to different real- 
izations of the stochastic velocity field. This condition is 
met sufficiently if the distance between two lines of sight is 
larger than one correlation length. A number of L/l x L/l 
lines of sight will then be adequate if the transverse ex- 
tent of the region is comparable to its depth. Figure 0] 
shows some composite spectra of 100 realizations, with 



R. Boger et al.: Effects of correlated turbulent velocity fields on the formation of maser lines 



5 




-10 o 10 



-10 o 10 



0.2 

o 
- — i 

X 

00 3 




-50 50 -50 50 

1.0 F ' ' "~ 






-10 10 




-10 10 

AvI Au D 




-50 50 




-10 10 




-10 10 



-10 10 



JL 




-10 10 




-50 50 -50 50 

0.9 



Az// Au D 

Fig. 4. Simulated spectra of an extended maser region. Each spectrum is composed of one hundred single realizations. 
In the upper ten panels we show the model calculations with a parameter set of l/L = 0.1, ot/uth — 10, tq = —110 
and in the lower ten panels the parameter set is l/L = 0.1, at/uth = ^0, tq — —400. The spectrum marked by an 
asterisk is used in Fig. |H1 



Table 1. For various correlation length the probability 
that the maximum intensity of an individual random spec- 
trum exceeds a certain value is given. For all cases the 
peak intensities are compared to the intensity expectation 
value in the line center as calculated from the parameter 
configuration given in Fig. [21 



Greater than 


L/l = 10 


L/l = 25 


L/l = 50 


0.1 ■ (Ja„=o) 


69.2% 


37.3% 


15.7% 


1 ■ (Ia»=o) 


28.0% 


3.5% 


0.13% 


10 • (Iau=o) 


6.5% 


0.06% 


0.0% 
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0.9% 


0.0% 


0.0% 



l/L = 0.1 and a t = 10 x v t h and er t = 50 x v t h respec- 
tively. It should be emphasized that these spectra are not 



selected by predefined criteria but are randomly chosen. 
It is conspicuous that they exhibit a distinct structure 
which is far from the spectrum of the intensity expec- 
tation value. They show pronounced maxima spreading 
partly over a frequency range of more than 2 x a%. In 
the standard analysis they would be attributed to several 
independent sources at different radial velocities. As men- 
tioned before, the spiky structure of the spectra is caused 
by the correlations in the velocity field. It should be noted, 
however, that - due to the exponential amplification - the 
intensity distribution does not reflect directly the veloc- 
ity distribution (see Fig. [2] in which the optical depth is 
proportional to the projected velocity distribution). 

Figure [S] shows a selected sample of observed OH and 
CH3OH maser spectra obtained by J. Caswell (private 
communication). A detailed description of most sources 
can be found in Caswell <jl1%?l l2UTTT|l together with 
additional maser spectra. Our simple model does clearly 
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Fig. 5. A sample of selected OH (left column) and 
CH 3 OH (right column) maser spectra observed by J. 
Caswell. All sources are associated with star forming re- 
gions or ultra compact Hn regions. Methanol masers are 
in general negligibly polarized which makes them partic- 
ularly suitable for a comparison with our synthetic spec- 
tra. Except of the two spectra of the source 309.921+0.479 
which are lefthanded polarized the polarization of all other 
OH sources is unknown 



not allow a detailed analysis of the observations. However, 
it is remarkable that it can reproduce the characteristic 
structure of the observed line profiles. Also the velocity 
scale of our model spectra is in agreement with the ob- 
servations, if for example, a typical cloud temperature of 
25 K is assumed. In this case the thermal velocity of the 
OH and CH 3 OH molecules is v th (OH) = 0.156 kms -1 
and Vt^CH^OH) = 0.114kms _1 , respectively. Due to the 
higher velocity dispersion, the spectra in the lower panels 
of Fig.0]show a richer structure and span a larger velocity 
range than the spectra in the top panels. Of course, also 
the amplification of radiation is affected by larger veloc- 
ity differences along a line of sight. Therefore, we had to 
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Fig. 6. Expectation values of the intensity calculated for 
an expanding plane parallel slap with a constant velocity 
gradient of dv sya /ds — 0, 5, 10, 15, 20, 25, 30 vth/L (in de- 
creasing order). All other parameters are the same as in 
Fig.H 



chose for both models very different opacity parameters 
to achieve the typical spiky structure of maser profiles. 

In this context, one may also ask how an additional 
systematic velocity field influences the line formation. For 
this purpose we considered an expanding plane parallel 
slab with a constant systematic velocity gradient super- 
posed on the stochastic velocity field. The general effect 
of this modification of our model can be viewed in Fig. El 
where the expectation value of the intensity for various 
values of the expansion velocity is shown (solid lines). As 
expected, this additional velocity component leads on av- 
erage to a weakening of the radiated emission and to a 
shift of the velocity distribution function. Additionally, a 
slight broadening of the profile of the intensity expecta- 
tion value is seen. Its width gives directly an estimation of 
the velocity range over which individual random spectra 
will be distributed. The shape of individual line profiles 
is very sensitive to changes of the velocity field. In Fig. [3 
we model again the radiation of a maser region by a su- 
perposition of L/l x L/l lines of sight. The turbulent ve- 
locity field is in all cases the same and only the strength 
of the linear expansion differs. From this example it is ob- 
vious that even small outflow velocities can significantly 
change the appearance of the resulting spectra. The ten- 
dency that with increasing systematic motions the emis- 
sion is reduced is also clearly seen. But the gradients of the 
systematic and the turbulent velocity field can also partly 
compensate each other, allowing a strong amplification of 
radiation along particular lines of sight. This behavior is 
seen in the sequence dv sys /ds = 0, 1, 2, 4 v t h/L and in the 
case dv sys /ds = 20 Vth/L. 

It is also instructive to simulate the spatial intensity 
distribution of composite spectra. As described before the 
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Fig. 7. As in Fig. 01 the spectra model the emission of an 
extended maser region with the parameter configuration: 
l/L = 0.1, ut/vth = 10, to = —110. For all calculations 
the same turbulent velocity field was used. Additionally 
we consider along the line of sights a systematic velocity 
component with a constant velocity gradient dv sys /ds. The 
values used in each calculation are shown in the upper left 
corner of the panels 
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Fig. 8. Spatial intensity distribution of the spectrum in 
FigQ] marked with an asterisk. Black circles indicate a 
maximum intensity of a single realization greater than ten 
times the expectation value in the line center. Grey and 
white circles show maximum intensities in the range be- 
tween 10 (I v = ) > Iv,max > (4=o) and (I V=Q ) > I v , max > 
0.1 (iy=o) respectively. The intensity in the blank boxes is 
below one tenth of the mean intensity in the line center 



emerging radiation from our model cloud is approximated 
by the superposition of L/l x L/l independent lines of 
sight corresponding to individual realizations of the ve- 
locity field. The general characteristics of the spatial in- 
tensity map of the cloud's surface can be simulated by 
distributing the maximum intensity found in each con- 
tributing spectrum randomly on a square grid of L/l x L/l 
cells. In Fig. [S] we consider as an example the spectrum 
marked with an asterisk in Fig. The black circles in- 
dicate maximum intensities greater than ten times the 
expectation value in the line center. Their number cor- 
responds to the number of distinct spectral components. 
The grey and white circles show maximum intensities 
in the range between 10 (I v =o) > Iv.m&x > (Iu=o) an d 
(Iv=o) > 4 ; max > 0.1 (ij/=o) respectively. The maximum 
intensity in the blank boxes is below one tenth of the ex- 
pectation value of the intensity in the line center. Due 
to the underlying stochastic velocity field, large spatial 
intensity variations are seen. Similar to observed high res- 
olution intensity maps the maser region consists of a few 
strong sources at different radial velocities (see, e.g., Rcid 
et al. I198()|l . while the radiation from most of the region 
is inconspicuous. 



4. Conclusion 

In the previous sections we have investigated the influence 
of a turbulent velocity field with finite correlation length 
on the appearance of a cosmic maser source. For this, we 
used a very simple stochastic model in which the only 
spatially varying quantity is the hydrodynamical velocity, 
while all other parameters, in particular the density and 
the pump efficiency are held constant. With this simple 
model, we could to a surprisingly high degree reproduce 
the observed characteristics of cosmic masers, in particu- 
lar those related to star forming regions. Our results are 
a strong indication that the details of the hydrodynamic 
velocity field inside a maser region are of great impor- 
tance for the interpretation of the observed properties of 
the source. Of course, reality is considerably more complex 
than our simple model. It is very likely that in star form- 
ing regions not only the velocity is a fluctuating quantity, 
but also the density and other physical parameters deter- 
mining the pump efficiency show strong spatial variations. 
However, our results convincingly show that the observa- 
tion of a maser spot, by itself does not necessarily imply 
that the physical conditions at the location of the spot are 
peculiar. The spot may as well be the result of the details 
of the velocity distribution along the line of sight. 

In view of the simplicity of our model, we have made 
no attempt to compare results in more detail with observa- 
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tional results of particular maser regions. Such attempts, 
however, have been performed by other authors investi- 
gating effects of turbulent motions on maser spectra. The 
approach closest to ours is that of Sobolev et al. (1998). 
They consider also a slab which is homogeneous except 
for the velocity field. The latter is characterized by its 
statistical properties. From these, individual realizations 
of the velocity field are constructed which then are used 
to calculate the optical depth and the emitted spectrum. 
Their approach differs from ours in the assumptions about 
the velocity field and in the numerical procedure. For the 
velocity field they assume a more general power spectrum 

ft(4) ~ «"W (16) 

where k mm has been introduced to avoid divergence for 
k — > 0. Physically, fc m i n is identified with the inverse of the 
linear size of the emitting volume. (In our case this identi- 
fication would correspond to the assumption l/L = 1.) To 
construct individual realizations of the velocity field ran- 
dom sampling is done in k-space (for details see Dubinski 
et al., 1995) and the velocity field is then obtained by a 
Fourier transformation. Sobolev et al. investigate in par- 
ticular models with a close to 17/6. The case a — 17/6 
corresponds to a Kolmogorov spectrum, which is the com- 
mon approximation for turbulence in an incompressible 
fluid. The general result obtained by Sobolev et al. is sim- 
ilar to ours: The spiky structure of maser spectra and the 
concentration of the emitted radiation in narrow bright 
spots, may be caused solely by the velocity field. - They 
then compare their results with detailed observations of 
the methanol masers in OMC-1. They find evidence that 
the power spectrum of the turbulence in OMC-1 is some- 
what steeper than an actual Kolmogorov spectrum and 
conclude that a > 3. (With a = 3 the spectrum (16) ap- 
proaches (7) in the limit k — ► oo) - We note in passing that 
for this more general velocity field there exists no general- 
ized radiative transfer equation as was derived by Traving 
and associates ( Gail et al. 1974, Traving 1975) for the spe- 
cial correlation function (3). I.e., if one is interested in the 
expectation value of the spectrum, one actually has to cal- 
culate the average of a sufficiently large number of spectra 
for individual realizations. In the case of maser emission 
this number may be very large, as shown in Fig. |3| 

A very different approach is followed by Gwinn 
(1994a, b). He tries to solve the inverse problem, i.e. he 
attempts to derive the characteristics of the velocity field 
directly from the very detailed observations available for 
the H2O maser region in W49N. The observations are 
interpreted by a model in which a strong stellar wind 
strikes ambient material and the arising shocks provide 
the excitation energy for the masers. More than 250 maser 
features were identified. Gwinn measures the distribution 
of the individual maser spots in coordinate and velocity 
space. From the derived power law correlation functions he 
concludes that the velocity field in the maser region is tur- 
bulent. This interpretation of the observations makes the 
(implicit) assumption that the velocities of the individual 



maser spots reflect directly the velocity field in the maser 
region. This assumption is appropriate if one considers the 
maser spots as physical entities like clumps. If one consid- 
ers models like ours or that of Sobolev et al. (1998), in 
which the maser spots are solely caused by the correla- 
tions in the velocity field, the situation is more complex 
(see Fig. 01. - A similar statistical analysis of high quality 
observational data for the H2O masers in the star form- 
ing regions Sgr B2(M), W49N, W51(MAIN), W51N, and 
W3(OH) has been performed by Strelnitski et al. (2002). 
They find that the two-dimensional distribution of maser 
spots shows a fractal structure and that the two-point ve- 
locity structure functions can be approximated by power 
laws with exponents close to Kolmogorov's values. From 
these findings they conclude that the velocity field in the 
these maser regions is highly turbulent. 
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